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Introduction 



The watershed transform is one of the major image segmentation tools [4], 
used in the community of mathematical morphology and beyond. If the 
watershed is a successful concept, there is another side of the coin: a 
number of definitions and algorithms coexist, claiming to construct a 
wartershed line or catchment basins, although they obviously are not 
equivalent. We have presented how the idea was conceptualized and 
implemented as algorithms or hardware solutions in a brief note :" The 
watershed concept and its use in segmentation : a brief history" 
(arXiv:1202.0216vl), which contains an extensive bibliography. See also 
|25| for an extensive review on the watershed concepts and construction 
modes. 
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Introduction 



The present work studies the topography of a relief defined on a node or 
edge weighted graph with morphological tools. The catchment basin of a 
minimum is defined as the sets of points it is possible to reach by a non 
descending path starting from this minimum. The watershed zone, i.e. the 
points linked with two distinct minima through a non ascending path, is 
more and more reduced as one considers steeper paths. 
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Outline 



• Reminders on weighted graphs 

9 Distances on node or edge weighted graphs 

• Adjunctions between nodes and edges 

• The flooding adjunction. 

• Invariants of the flooding and closing adjunction 

• Flooding graph 
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Outline 



a Paths of steepest descent and k-steep graphs 

• The scissor operator, minimum spanning forests and watershed 
partitions 

• Lexicographic distances and SKIZ 

• The waterfall hierarchy. 

• Emergence and role of the minimum spanning tree 

• Discussion and conclusion 
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Reminder on adjunctions 



The following section contains a reminder on adjunctions linking erosions 
and dilations by pairs and from which openings and closings are derived 
ESI.fia.fll] 
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Preamble on adjunctions 



Let T a complete totally ordered lattice, and T>,£ be arbitrary sets 

: the smallest element and Q : the largest element of T 
Fun(X',T) : the image defined on the support V with value in T 
Fun(£^,T) : the image defined on the support £ with value in T 

Let f be a function of Fun(I',T) and g' be a function of fun{£,T) 

01 : ¥un{V,T) ^ Fun(£:,r) and 

j6 : Fun(£^,T) — ?■ Fun(X',T) be two operators 



Definition 



OL and /3 form an adjunction if and only if : 

for any f in Fun(I?,T) and g in Fun(£^,T) ; ocf < g <^ f < 
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Erosions and dilations 



Theorem 



If (a, fi) form an adjunction, then <x is a dilation (it commutes with the 
supremum of functions in Fun(X',T) ) 

and j6 is an erosion (it commutes with the infimum of functions in 

Fun(f,r) ; 

Proof: Let f be a function of Fun(X',T) and {g)i functions of Fun(£^,T) 
Then 

f < ^Agi ^af < Agi O V/ : af < g; ^^i -.f < ^g; ^ f < A^g; 

/ / / 

Reading these equivalences from left to right and replacing f by ^f\gi 

implies that ^Ag; < AjSg"/. 

Reading these equivalences from right to left and replacing f by Ajigi 

implies that A^g; < ^Agi 

i / 
establishing that ^/\gi = A^g-,. 
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Erosions and dilations (2) 



Similarly we show that ^ is a dilation, i.e. commutes with the union. 
Remark: Calling a constant function equal to 0, we have for any f : 
< f>f implying aO < f. This relation is true for all f, indicating that 
uO = 0. 

Similarly, we show that j5Cl = Q, where Q is the constant function equal 
to Q. 
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Increasing operators 



Lemma 



DL and /3 are increasing operators. 



Proof. 



\ff<g then li{fAg) = fi{n = fiin A ^{g) < ^{g) 



□ 
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From one operator to the other 



From the relation af < g ^ f < f)g one derives expressions of one 
operator in terms of the other one. 

If DC is known, j6 may be expressed = \/ {f \ af < g} . 

Inversely if f> is known, then ocf = f\\^g \ f < ^g} 
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Opening and closing 



Lemma 



The operator is a closing : increasing, extensive and idempotent. 
Similarly the operator ocji is an opening : increasing, anti-extensive and 
idempotent. 

■ 

Proof: Openings and closings, obtained by composition of increasing 
operators, are increasing. 

By adjunction we obtain af < af ^ f < jiaf showing that the closing is 
extensive and ocjif < f <J= /3f < /3f showing that the opening is 
anti-extensive. 

In particular, applying an opening to af yields ocf > ot^otf . 
On the other hand, a being increasing, and the closing extensive : 
f < fiaf =^ af < a.fia.f showing that Kf = afiaf 
Applying j6 on both sides yields fiaf = ^oc^ocf 
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The family of invariants of an opening or a closing 



We call lnv(7) and lnv((p) the family of invariants of 7 and cp. 
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The family of invariants of an opening 



Lemma 



The family of invariants of an opening is closed by union. And the family 
of closings is closed by intersection. 

Proof: Let us prove it for openings ; the result for closings being obtained 
by duality. Suppose that gi and g2 are invariant by the opening 7 : 
7{gi) = gi and 7(5-2) = g2- 

Then j{giV g2) < giV g2 = j{gi) V 7(5-2) by antiextensivity 
And 7(51 Vg-2) > l{gi) V7(g-2) as 7 is increasing. 
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Invariants of an opening or a closing 



9 ocf ^ lnv(7) as ocf = afiaf = (a/5)af 
» e Inv(^) as = ^a^f = {^a)^f 
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The opening as pseudo-inverse operator of the erosion 



* The erosion has no inverse as al^g < g. 

* The opening is the pseudo inverse of the erosion : alig is the smallest 
set having the sanne erosion as g. 

Proof : Suppose that f verifies fif = fig which implies Oi^f = af^g. if 
f < we have afif < f < oc^g showing that f = a^g 

* if If G lnv(7), we have ^g = a^g = g : on lnv(7), oc = fi^^ 
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The closing as pseudo-inverse operator of the dilation 



* The dilation has no inverse as jiag > g. 

* The closing is the pseudo inverse of the dilation : jiag is the largest set 
having the same dilation as g. 



* if If G Inv(^), we have cpg = fiag = g : on \nv{(p), j6 = a 



-1 
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Reminder on graphs 
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Graphs : General definitions 



A non oriented graph G = [N , E] : N = vertices or nodes ; E = edges ; 
an edge u ^ E = pair of vertices (see [I].|10]) 

A chain of length n is a sequence of n edges L = {tvi, U2, . . . , Un}, such 
that each edge u, of the sequence (2 < /' < n — 1) shares one extremity 
with the edge u/_i (t7;_i u-,), and the other extremity with [v/+i 

^ Ui). 

A path between two nodes x and y is a sequence of nodes 

(ni = X, n2, n/t = y) such that two successive nodes n; and n/_|_i are 

linked by an edge. 

A cycle is a chain or a path whose extremities coincide. 



A cocycle is the set of all edges with one extremity in a subset Y and the 
other in the complementary set Y. 

I II i| I mil 
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Graphs : partial graphs and subgraphs 



The subgraph spanning a set A C N : Ga = [A, Ea], where Ea are the 
edges linking two nodes of A. 

The partial graph associated to the edges E' C E '\s G' = [N, E'] 
For contracting an edge (ij) in a graph G, one suppresses this edge u 
and its two extremities are merged into a unique node k. All edges 
incident to / or to J become edges incident to the new node k. 
Contraction : 

If /-/ is a subgraph of G, The operator k : {G , H) ^ G' contracts all edges 
of AV in G 



cr 
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Graphs : Connectivity 



A connected graph is a graph where each pair of nodes is connected by a 
path 

A tree is a connected graph without cycle. 
A spanning tree is a tree containing all nodes. 
A forest is a collection of trees. 

Labelling a graph = extracting the maximal connected subgraphs 
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Weighted graphs 



In a graph [N, E], N represents the nodes, E represents the edges and T 
is the set of weights of edges or nodes. (0 being the lowest weight and Q 
the largest) 

Edges and nodes may be weighted : ejj is the weight of the edge (/',_/') and 

Pi the weight of the node /'. 

The following weight distributions are possible: 

( — , o) : no weights 

( — , n) : weights on the nodes 

(e, o) : weights on the edges 

(e, n) : weights on edges and nodes 

Various types of graphs are met in image processing. 
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Which types of graphs ? 



In "pixel graphs" the nodes are the pixels and the edges connect 
neighboring pixels. The weights of the pixels are their value and the 
weights of edges may be for instance a gradient value computed between 
their extremities. 

O^O^O 

\//\^\^\^\/^\// 

o^o^o 
• — o^o^o 

\/\/\//\/\/ 
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Which types of graphs ? 



We will work with "neighborhood graphs" where the nodes are the 
catchment basins and the edges connect neighboring bassins. The edges 
are weighted by a dissimilarity measure between adjacent catchment 
basins; the simplest being the altitude of the path-point between two 
basins. 
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The gabriel graph 



Gabriel and Delaunay graphs permit to define neighborhood relations 
between the sets and nodes of a population. 




Two nodes x and y are linked by an edge if there exists no other node in 
the disk of diameter (xy). Node weights will express features of the nodes, 
edge weights relationships between adjacent nodes. 
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Gabriel graph modeling a lymphnode 



The nuclei of a lymphnode are segmented and their Gabriel graph 
constructed. 



• 1.1 : histological section. 

• 1.2 : marking the nuclei. 

• 2.1 : catchment basins of 
the nuclei on the original 
image. 

• 2.2 : node weighted 
graph: the node weights 
represent the nucleus 
type. 
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Minimum spanning trees and forests in an edge weighted 
graph 



A minimum spanning tree is a spanning tree for which the sum of the 
edges is minimal. 

A spanning forest is a collection of trees spanning all nodes. 

In a minimum spanning forest, the sum of the edges is minimal. 

The minimum minimorum spanning forest (MSF) has only nodes. Wih an 

additional constraint, one gets particular forests: 

• a MSF with a fixed number of trees, useful in hierarchical 
segmentation [T^ 

• a MSF where each tree is rooted in predefined nodes, useful in marker 
based segmentation. 
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Flat zones and regional minima on edge weighted graphs 



A subgraph G' of an edge weighted graph G is a flat zone, if any two 

nodes of G' are connected by a chain of uniform altitude. 

A subgraph G' of a graph G is a regional minimum if G' is a flat zone and 

all edges in its cocycle have a higher altitude 

fie : G ^ jigG : extracts all regional minima of the graph G 



Flat zones and regional minima on node weighted graphs 



A subgraph C of a node weighted graph G is a flat zone, if any two nodes 
of G' are connected by a path where all nodes have the same altitude. 
A subgraph G' of a graph G is a regional minimum if G' is a flat zone and 
all neighboring nodes have a higher altitude 
Hr,: G ^ finG : extracts all regional minima of the graph G 



Contracting or expanding graphs 



Contraction: 

For contracting an edge (ij) in a graph G, one suppresses this edge u and 
its two extremities are merged into a unique node k. All edges incident to 
/ or to j become edges incident to the new node k and keep their weights 

If /-/ is a subgraph of G, The operator k : [G , H) ^ G' contracts all edges 
of AV in G 

Expansion: 

The operator ^ : ( — , n) — )■ ^ G creates for each isolated regional 
minimum / a dummy node with the same weight, linked by an edge with /. 
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Extracting partial graphs 



The following operators suppress a subset of the edges in a graph: 
X : (— , o) — )■ keeps for each node only one adjacent edges. 
J, : (e, o) — > J, G keeps for each node only its lowest adjacent edges. 
J| : ( — , n) ^ J| G keeps only the edges linking a node with its lowest 
adjacent nodes. 

These operators may be concatenated: 

^ keeps for each node only one lowest adjacent edge. 
X \ keeps for each node only one edge linking it with its lowest adjacent 
nodes. 
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Distances on a graph 

Case of edge weighed graphs 
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Constructing distances on an edge weighted graph. 



Distances on an edge weighted graph have chains as support : 

1) Definition of the weight of a chain, as a measure derived from the edge 
weights of the chain elements (example : sum, maximum, etc.) 

2) Comparison of two chains by their weight. The chain with the smallest 
weight is called the shortest. 

The distance c/(x,y) between two nodes x and y of a graph is oo if there 
is no chain linking these two nodes and equal to the weight of the shortest 
chain if such a chain exists. 

Given three nodes (x,y,z) the concatenation of the shortest chain Uxy 
between x and y and the shortest chain Uyz between y and z is a chain 
TTxz between x and z, whose weight is smaller or equal to the weight of 
the shortest chain between x and z. To each distance corresponds a 
particular triangular inequality : c/(x, z) < weight{ 7Zxy\> TCyz) where 
nxy\> Tiyz represents the concatenation of both chains. 
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Distance on an edge weighted graph based a the length of 
the shortest chain 



Length of a chain: The length of a chain between two nodes x and y is 
defined as the sum of the weights of its edges. 
Distance: The distance d{x,y) between two nodes x and y is the 
minimal length of all chains between x and y. If there is no chain between 
them, the distance is equal to oo. 

Triangular inequality : For (x,y,z) : d{x,z) < d{x,y) + d{y,z) 
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Distance on a graph based on the maximal edge weight 
along the chain 



The weights are assigned to the edges, and represent their altitudes. 
Altitude of a chain: The altitude of a chain is equal to the highest 
weight of the edges along the chain. 

Flooding distance between two nodes: The flooding distance 
fldist(x,y) between nodes x and y is equal to the minimal altitude of all 
chains between x and y. During a flooding process, in which a source is 
placed at location x, the flood would proceed along this chain of minimal 
highest altitude to reach the pixel y. If there is no chain between them, 
the level distance is equal to oo. 

Triangular inequality : For (x,y,z) : c/(x,z) < cl{x, y) V d{y, z) : 
ultrametric inequality 
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The flooding distance is an ultrametric distance 



An ultrametric distance verifies 

* reflexivity : c/(x,x) = 

* symmetry: d{x,y) = d{y,x) 

* ultrametric inequality: for all x, y, z : c/(x, y) < max{d{x, z), d{z, y)} : 
the lowest lake containing both x and y is lower or equal than the lowest 
lake containing x, y and z. 



Distances on a graph : sum and maximum of the edge 
weights 




The shortest chain (sum of weights of the edges) between x and y is a red 
line and has a length of 4. 

The lowest chain (maximal weight of the edges) between x and y is a red 
line and a maximal weight of 2. A flooding between x and y would follow 
this chain. 

For this particular example, the shortest and lowest paths are identical. 
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The lexicographic distance or the cumulative effort for 
passing the highest edges 



Toughness of a chain: We call toughness of a chain the decreasing list of 
altitudes of the highest edges met along this chain. 
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In fig. A, along the bold chain between x and y, the highest edge rises at 
the altitude 3. After crossing it, the highest edge on the remaining chain 
rises at 2. After crossing it, one is at destination. Hence the toughness of 
the chain from x to y is [3,2]. Distances are compared in a lexicographic 
order. The shortest chain is the red chain in fig.B, it has two edges with 
weight 2 and its toughness is [2,2]. pi t|| | |||||| | || | || 
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The lexicographic distance : a more formal definition 



The lexicographic length A{A) of a chain A = ei2^3---Sn-in is 
constructed as follows. Following the chain from the origin towards its end, 
one records the highest valuation of the chain, let it be Ai, then again the 
highest valuation A2 on the remaining part of the chain and so on until the 
end is reached. One gets like that a series of decreasing values: 
Ai > A2 >....> An- 

The lexicographic distance between two nodes x and y will be equal to the 
shortest lexicographic length of all chains between these two nodes and is 
written lexdist(x, y). 

Remark: The lexicographic length of a never increasing path (NAP), is 
simply equal to the series of weights of its edges, as it is the highest edge 
along the lowest path between x and y. Later we introduce shortest path 
algorithms whose geodesies are NAPs. 
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The lexicographic distance of depth k 



The lexicographic distance of depth k between two nodes x and y is 
written lexdist/((x, y) and obtained by retaining only the k first edges. in 
lexdist(x,y). 



Remark: lexdisti(x, y) is the same as the ultrametric flooding distance as 
it is the highest edge on the lowest path between x and y. 
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Distances on a graph 

Case of node weighed graphs 
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Constructing distances on a node weighted graph. 



Distances on a node weighted graph have paths as support : 

1) Definition of the "length" of a path, as a measure derived from the 
node weights of the path elements (example : sum, maximum, etc.) 

2) Comparison of two paths by their length. The path with the smallest 
length is called the shortest. 

The distance c/(x,y) between two nodes x and y of a graph is oo if there 
is no path linking these two nodes and equal to the length of the shortest 
path if such a path exists. 

Given three nodes (x,y,z) the concatenation of the shortest path Uxy 
between x and y and the shortest path Tiyz between y and z is a path tCxz 
between x and z, whose length is smaller or equal to the length of the 
shortest path between x and z. To each distance corresponds a particular 
triangular inequality : d{x,z) < length{ TTxy> Jiyz) where TCxy\> TCyz 
represents the concatenation of both paths. 
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Distance on a node weighted graph based on the maximal 
node weight along the path 



The weights are assigned to the nodes, and represent their altitudes. 
Altitude of a path: The altitude of a path is equal to the highest weight 
of the nodes along the path. 

Flooding distance between two nodes: The flooding distance 
fldist(x,y) between nodes x and y is equal to the minimal altitude of all 
paths between x and y. 



Distance on a node weighted graph based on the cost for 
travelling along the cheapest path 



The weights are assigned to the nodes and not to the edges. Each node 

may be considered as a town where a toll has to be paid. 

Cost of a path: The cost of a path is equal to the sum of the tolls to be 

paid in all towns encountered along the path (including or not one or both 

ends). 

Cost between two nodes: The cost for reaching node y from node x is 
equal to the minimal cost of all paths between x and y. We write 
tolldist(x, y). If there is no path between them, the cost is equal to oo. 



I 



ustration of the cheapest path 













































In this figure, the cheapest chain between x and y is in red and the total 
toll to pay is 1 + 1 + 2 + 2 = 6 



I' 1 Illllllllll II I lllll 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigi p012 February 29 45 / 201 



Two adjunctions between edges and nodes on weighted 

graphs 
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Two adjunctions between edges and nodes 



Definition 



We define two operators between edges and nodes : 

- an erosion [eeni]/,- = A nj and its adjunct dilation 

[Snee]i= y eik 

{k neighbors of /) 

- a dilation [Senn] - = n-, V nj and its adjunct erosion 
[£nee],- = ^eik 

{k neighbors of /) 



Lemma 



The operators we defined are pairwise adjunct or dual operators: 
- £ne 3nd Sen 3re adjunct operators 

£en 3nd 5ne sre adjunct operators 

Ene 3nd Sne are dual operators 

ten and Sen are dual operators 



I ji ij I II II 11 1 II t). 
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Let us prove that 3en and Sne are adjunct operators 

If G = [e, n] and G = [e, n] are two graphs with the same nodes and 
edges, but with different valuations on the edges and the nodes then 

Senn <e^ yij : n; V nj < e;j <^ yij : n, < e;j <^ yij : n, < 

ASy = [enee]i O /I < EneS 

{j neighbors of /) 

which establishes that Sen and £en are adjunct operators. 

Let us prove that £en and Sen are dual operators 

[£en{-n)];j = -n, A -nj = -{n, V nj) = - [denn]ij 
Hence Senn = — [Sen ( — ")] 
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The flooding adjunction 



The dilation [(5enn]y = n, V nj and its adjunct erosion 

[Ene^]/ = A s//< have a particular meaning in terms of flooding. 

{k neighbors of /) 

If rij and nj represent the altitudes of the nodes /' and j, the lowest flood 
covering / and J has the altitude [Senn]ij = n, V nj 

If /' represents a catchment basin, e//( the altitude of the pass points with 
the neighboring basin k, then the highest level of flooding without 
overflow through an adjacent edge is [£nee],- = A s/Zc 

{k neighbors of /') 
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Waterfall flooding on graphs : an erosion 



The waterfall flooding is completely specified if one knows the level of 
flood in each catchment basin 



• The waterfall flooding fills 
each catchment basin up 
to its lowest pass point. 

• In terms of graphs : the 
flooding in a node is 
equal to the weight of its 
lowest adjacent edge. 

It is an erosion between 
node weights and edge 
weights : 

[enee]i= Ae/fc 

(k neighbors of /) 




a bed e f g h I | 

Erosion from the edges to the nodes 



llllllllll II ill! 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigl 



012 February 29 



Waterfall flooding on graphs : an erosion 



The waterfall flooding is completely specified if one knows the level of 
flood in each catchment basin 



The waterfall flooding fills 
each catchment basin up 
to its lowest pass point 
In terms of graphs : the 
flooding in a node is 
equal to the weight of its 
lowest adjacent edge 

It is an erosion between 
node weights and edge 
weights : 

[enee]i= Ae/fc 

(k neighbors of /) 




I 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigl 



012 February 29 



Erosion between edges and nodes 




Dilation between nodes and edges 



If n, and nj represent the altitudes of the nodes / and j, the lowest flood 
covering / and j has the altitude [Senn] - = n, V nj 




dilation from nodes to edges 
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Erosion between edges and edges / between nodes and 
nodes 



By concatenation of operators between edges and nodes we obtain : 

• an adjunction between nodes and nodes : {Sne^en: ^ne^en) = (Sn.^n) 

• an adjunction between edges and edges ; {een^ne, ^en^ne) = (Ce.^e) 
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Opening and closing 
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Opening and closing 



As £ne and Sen are adjunct operators, the operator (pn = Cne^en is a closing 
on n and 7e = ^en£ne is an opening on e 

Similarly the operator fe = £enSne is a closing on e and 7n = Sne£en is an 
opening on n 

In the sequel, the flooding adjunction will play a key role, in particular the 
associated opening 7e and closing cpn- 
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The opening je 
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I 



ustration in one dimension of the opening 7^ 




Erosion from the edges to the nodes 




Dilation from the nodes to the edges 



Sen^ne ■ on the left, the result of the erosion, filling each basin to its 



lowest pass point. On the right the subsequent dilation. Some pass points 
have a reduced altitude : the amount of reduction is indicated in red. 
These passpoints are those which are not the lowest pass points of a 
catchment basin. 
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ustration on an edge weighted tree of the opening je 




7e = SenSne ■ From left to right: 1) an edge weighted graph, in the centre, 
2) the result of the erosion e^e, 3) the subsequent dilation produces an 
opening. The edges in red are those whose weight has been reduced by the 
opening. These edges are not the lowest edges cfpj| n; 'j i 'l' l 'ai [| [ d[|L | | [| t|n i| iL i tjb| [Yi 
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The invariants of the opening je 



Two possibilities exist for an edge with a weight A : 

• the edge (ij) has lower neighboring edges at each extremity. Hence 

ene{i) < A and e^eO') < A ; hence 3enene{i,J) = een(') V£en(j) < A 

• the edge (/,_/') is the lowest edge of the extremity /'. Then ene(') = A 

and SneU) < A ; hence SenSneiiJ) = £en{i) V SenU) = A 

Conclusion: the edges invariant by the opening 7e are the edges which are 
the lowest edge of one of their extremities. All edges with lower adjacent 
edges at their extremities have their weight lowered by the opening 7e 
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Extracting from an edge weighted graph a partial graph, 
invariant by the opening je 



The relation Sne = £ne{^en£ne) shows that all edges which are not 
invariant by the opening 7e = Sen^ne play no role in the erosion. As a 
matter of fact, if the opening lowers the valuation of an edge and if 
the subsequent erosion Cne is not modified, it means that the presence of 
this edge with its weight plays no role in the erosion e^e- 
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The waterfall graph, partial graph of the lowest adjacent 
edges 



Suppressing in an arbitrary graph g all edges which are not invariant by 
the opening 7e produces a graph g' , invariant for the opening 7e and 
called waterfall graph. The operator keeping for each node only its lowest 
adjacent edges is written J, : (e, o) — > J, G 
Properties: 

• J, G spans all the nodes (each node has at least one lowest 
neighboring edge ; such edges are invariant by 7e) 

• And £ne(G) = £f)e(i <j) {tne assigns to each node the weight of its 
lowest adjacent edge, which is the same in G as in J, G) 
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Illustration on a planar edge weighted graph of the opening 

7e- 




erosion opening waterfall graph = edges 

invariant by opening 
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Erosion from edges to nodes on the graph 4- G 



1 1 




A: Initial graph G 

B: Erosion £ne from the edges to the nodes on G 
C: The graph J, G 

D: Erosion £ne from the edges to the nodes on J, i^j: " ( '"^) ||-|- || fli|je|(-l|| *^ | ) ]r[n 
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Erosion from edges to nodes on the graph 4- G 



Another illustration of £ne[G) = £ne{i G) 




erosion on the RAG erosion on the waterfall graph 
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The invariants of the opening je 



I' 1 II I mill I II I II II I 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigi p012 February 29 66 / 201 



The invariants of the opening je 



If (^1.^) and {e2,o) are two edge weight distributions which are invariant 
for 7e, then (ei V 62,0) also is invariant for 7e. 

An arbitrary graph may be transformed into a flooding graph: 

* For an arbitrary edge weight distribution (e,o) : (7ee,o) G lnv(7e) 

* For an arbitrary node weight distribution ( — , n) : {Senn, n) G lnv(7e) as 
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The invariants of the opening je 



For a connected graph g = (e,o) G lnv(7e): 
9 any partial spanning graph belongs to inv(7e) 
• any subgraph belongs to lnv(7e) 

In particular J, : G = (e, o) — > J, G containing for each node only its 
lowest adjacent edges 

and X i G keeping only one lowest adjacent edge for each node, 
both belong to lnv(7e) 
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The regional minima of the opening jei 



Theorem 



If G = (e, o) G lnv(7e) and m = (e, o) is the subgraph of its regional 
minima, then Sne^ = {—,£nes) is the subgraph of the regional minima of 
the graph SneG = (-, Enee) 

Proof: A regional minimum m/< of the graph G = (e,o) G lnv(7e) is a 
plateau of edges with altitude A, with all external edges having a weight 
> A. If a node / belongs to this regional minimum, its adjacent edges have 
a weight > A but it has at least one neighboring edge with weight A : 
hence enee{i) = A. Consider now an edge (s, t) outside the regional 
minimum, with the node s inside and the node t outside the minimum. 
Then est > A. As G = (e, o) G lnv(7e), the edge (s, t) is then one of the 
lowest edges of the nodes t : thus £ne{t) = est > ^- This shows that the 
nodes spanned by the regional minimum m^ form a regional minimum of 
the graph £neG = { — ,£nee) 
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Inverse of £ne on the invariants of the opening je 



On lnv(7e) : 5en^ne = Identity showing that on lnv(7e) := = £ 
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The closing (pn 
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The closing (pn 



The closing fn is obtained by a dilation Sen of the node weights followed 
by an erosion £ne- One remarks on the following figure that the node 
weights remain the same, except the isolated regional minima, which take 
the weight of their lowest neighboring node. We give the proof in the next 
slide. 
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The invariants of the closing cpn 



Two possibilities exist for a node / with a weight A : 

• the node / is an isolated regional mininnum. Then Sen assigns to all 
edges adjacent to / a weight bigger than A. The subsequent erosion 
£ne assigns to / the smallest of these weights. 

• the node /' has a neighbor J with a weight }i < A. Then Sen assigns to 
the edge (/',_/) the weight A; whatever the weight of the other 
adjacent edges. The subsequent erosion ene assigns to / the smallest 
of these weights, that is A. 

The closing cpn replaces each isolated node constituting a regional 
minimum by its lowest neighboring node and leaves all other nodes 
unchanged. 
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The invariants of the closing (pn 
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The invariants of the closing cpn 



A node weighted graph is invariant for the closing cpn iff it does not 
contain isolated regional minima. 

For an arbitrary node weight distribution g = ( — , n) : —o g creates for 
each isolated regional minimum /' a dummy node with the same weight 
linked by an edge with /'. Hence —og^ lnv((pn) 



The invariants of the closing cpn 



If ( — , ni) and ( — , ni) are two node weight distributions which are 
invariant for cpn, then ( — , ni A /i2) also is invariant for fn- 
For an arbitrary node weight distribution ( — , n) : ( — , cpn^) G Inv(^n) 
For an arbitrary edge weight distribution (e,o) : (e.enes) G lnv((pn) as 
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The invariants of the closing cpn 



For a graph g = (e,o) G \nv{<pn): any partial or subgraph which does not 
create an isolated regional minimum also belongs to Inv(^n) 

Each node / belonging to a regional minimum of g is linked with at least 
another node j in this minimum through an edge with the same weight. 
This edge is one of the lowest adjacent edges of / and links / with one of 
its lowest neighboring nodes. 

For this reason after pruning, the graphs J, g and J| g still belong to 
lnv((pn). 
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The regional minima of the closing cpn 



Theorem 



If G — ( — , n) E Inv(^n) and m = ( — , n) is the subgraph of its regional 

minima, then Senm = {Senn,o) is the subgraph of the regional minima of 

the graph SenG = {3enn, o) 

— ^ 

Proof: A regional minimum m; of a graph G = { — ,n) G lnv((pn) is a 

plateau of pixels with altitude A, containing at least two nodes (there are 

no isolated regional minima in Inv(^n)). All internal edges of the plateau 

get the valuation A by Senn. If an edge (/,_/') has the extremity / in the 

minimum and the extremity j outside, then Senn{i,j) > A. Hence, for the 

graph {Senn, o), the edges spanning the nodes of m,- form a regional 

minimum. 
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Inverse of Sen on the invariants of the closing (pn 



On Inv(^n) : £ne^en = Identity showing that on lnv((pn) : £ne = ^, 
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The flooding graphs 
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The flooding graph 



Definition 



An edge and node weighted spanning graph G = [N, E] is a flooding 
graph iff its weight distribution (n, e) verify the relations: 

- Senn = e 

- EneS = n 



Corollary 



For a flooding graph weight distribution (n, e) : 

- n E lnv(g?n) 

- e G lnv(7e) 



Proof. 



e = 5enn = Sen^ne^ = le^ and n = £nee = Sne^enn = (fn^ □ j 
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Properties of the flooding graph 



As G is invariant by 7e, all its edges are the lowest edge of one of their 
extremities. 

As G is invariant by tpn, it has no isolated regional minimum. 
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The lowest adjacent edges of each node in a flooding graph 



In a flooding graph, £nee = n, hence all edges adjacent to a node have 
weights which are higher or equal than this node and at least one of them 
has the same weight. 

On the other hand, each node / has at least one neighbor j which is lower 
or equal (otherwise it would be an isolated regional minimum). The weight 
e,y verifies e,y = 5enn{i,j) = n,- V rij = n;. This shows that the edges linking 
a node with lower or equal nodes have the same weight than this node. 

In particular, the edges linking a node / to its lowest neighboring nodes 
belong to the lowest adjacent edges of this node and have the same 
weight: J| G C ^ G and G =J|^ G = J| G 
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Constructing flooding graphs. 



If G = (e,o) G lnv(7e) then (e.g^gs) is a flooding graph (since 
If G = ( — , n) G Inv(^n) then {3enn, n) is a flooding graph (since 

n = Cpnn = Sne^enn = Enee) 
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Deriving flooding graphs from ordinary graphs 



If G = (e,o) is an arbitrary edge weighted graph, \. (e,o) 6 lnv(7e), as 
the edges lowered by 7e, have been suppressed, the other edges keeping 
their weights. Recall that = £ne i e. The derived flooding graph 
simply is (| e, Cne I e) 

If G = (— , n) is an arbitrary node weighted graph, ( — , —on)^ Inv(^n), as 
isolated regional minima, if any, have been duplicated. The derived 
flooding graph simply is {Sen —°n,—o n). 
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Partial graph of a flooding graph 



Suppressing edges in a flooding graph, but leaving at least one lower 
neighboring edge for each node (like that, no isolated regional minima are 
created) produces a partial graph which also is a flooding graph, with the 
same distribution of weights on the nodes and on the remaining edges. 
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Regional minima of flooding graphs 
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Regional minima of a flooding graph 



We proved earlier these theorems: 

• If G = (e, o) G lnv(7e) and m = (e, o) is the subgraph of its regional 
minima, then Snem = (— , Enes) is the subgraph of the regional 
minima of the graph SneG = { — ,£nee). 

• If G = ( — , n) G lnv((pn) and m = ( — , n) is the subgraph of its 
regional minima, then Senin = {5enn,o) is the subgraph of the 
regional minima of the graph SenG = {Senn,<>). 

As in a flooding graph n = £neS and e = Senn we derive: 



Theorem 



If G is a flooding graph with the weight distribution (e, n), then the node 
weighted graph ( — , n) and the edge weighted graph (e, o) have the same 
regional minima subgraph 
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The regional minima on the edge or node graph within the, 
flooding graph 




regional minima on the edges regional minima on the nodes 
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Labeling the regional minima on the edge or node graph 
within the flooding graph 



As the minima are identical on the nodes or the edges of a flooding graph, 
it is possible to assign the same labels to nodes or to edges. 




labeling the regional minima 
on the nodes and/or the edges 
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Paths of steepest descent and catchment basins 
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The catchment basins of the minima 



Lemma 



From each node outside a regional minimum starts a never ascending path 

to a regional minimum in a flooding graph. 

— M 

Proof: Any node / outside a regional minimum has a lower neighboring 

node, if it does not belong to a plateau. Otherwise it belongs to a plateau, 

containing somewhere a node J with a lower neighboring node outside, as 

the plateau is not a regional minimum. The plateau being connected, there 

exists a path of constant altitude in the plateau between /' and j. Following 

this path, it is possible, starting at node / to reach the lower node k. This 

shows that for each node there exists a never ascending path to a lower 

neighboring node. Taking this new node as starting node, a still lower 

node may be reached. The process may be repeated until a regional 

minimum is reached. 

Thanks to this lemma, it is possible to define the catchment basins of the 
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The catchment basins of the minima 




The catchment basin of a minimum m is the set of all nodes from which 
starts a never ascending path towards m. 

As in a flooding graph, each node and its lower neighboring edges have the 
same weight, each node in such a never ascending path is followed by an 
edge with the same weight except the last one belonging to a regional 
minimum. For this reason the catchment basins based on the node weights 
or on the edge weights are identical. 

The watershed zones are the nodes belonging to more than one 
catchment basin. 

The restricted catchment basins are the nodes which belong to only 
one catchment basin : it is the difference between the catchment basin of 
a minimum m and the union of catchment basins of all other minima. 
From a node in a restricted catchment basin, there exists a unique non 
ascending path towards a unique regional minimum 




M-flooding graphs 



The catchment basins rely entirely on the non ascending paths of the 
graph reaching a regional minimum. The altitude of the regional minimum 
has no importance. If we consider only the end points of such a path, the 
fact that the minimum is an isolated node or not has no importance either. 
For this reason, we may relax the definition of the flooding graphs for 
which Senn = s and CneS = n for all edges and nodes. 
Consider a flooding graph where all non regional minima nodes have a 
positive weight. Assigning to all minima a weight does not invalidate the 
relation Senn = e. The relation Ene^ = n also remains true, except for 
isolated regional minima. We call M-flooding graphs, the node and edge 
weighted graphs verifying Senn = e everywhere and tneS = n for all nodes 
which are not isolated regional minima. 

In what follows we consider NAP which end with a node in a regional 
minimum. Whether this minimum is isolated or not has no importance. 
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Extending the restricted cathment basins and reducing the 
watershed zone : steep, steeper, steepest flooding graphs 
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Catchment basins and segmentation. 



The watershed transform is mainly used for segmentation, with the aim to 
create a partition representing precisely the extension of each object. Large 
watershed zones are ambiguous as they separate restricted catchment 
basins without precise localisation of the contour separating them. 
For obtaining precise segmentations, it is important to reduce these zones 
and even suppress them completely if possible. Reducing the number of 
never ascending paths from each node to a regional minimum would help 
constraining the construction of the watershed partition. Ideally, if only 
one such path remains for each node outside the regional minima, the 
solution would be unique, the restricted catchment occupy the whole 
space and the watershed zones be empty. 

In this section we show how to reduce the number of paths and keep only 
those which have some degree of steepness. Increasing the steepness 
reduces the number of paths. 
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Flooding tracks 



In a flooding graph each edge is the lowest edge of one of its extremities 
and has the same weigh: we call such a pair made of a node and adjacent 
edge with the same weight flooding pair. 

Two couples (/', ij) and (JJk) are chained if the node in the second couple 
is an extremity of the edge of the first couple. 

A flooding track is a list of chained couples of never increasing weight. 
The lexicographic weight of a flooding track is the list of never 
increasing weights of its pairs. 

Two flooding tracks may then be compared by comparing their 
lexicographic weights using the lexicographic order relation. 
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Pruning the flooding graph to get steeper paths. 



We define a pruning operator J,'' operating on a flooding graph G. The 
pruning J,'* considers each node outside the regional minima and 
suppresses its adjacent edges, if they are not the highest edge of a flooding 
track verifying: 

• their lexicographic weight is minimal. 

• their length is k. It may be shorter if its last couple belongs to a 
regional minimum 

After pruning, each node outside the regional minima is the origin of one 
or several k-steep flooding tracks (remark that the pruning only suppresses 
the highest edge of the track). If there are several of them, they have the 
same weights. 

We call them the k-steep adjacent paths to the node /'. We say that the 
graph J,'' G has a k-steepness or is k-steep. 
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Nested k-steep graphs 



As the steepness degree increases, the pruning becomes more and more 
severe, producing a decreasing series of partial graphs, hence for 
/c> / : G c;' G. Furthermore G =i^i^ G =1''^/ c 

The pruning J,^ does nothing as each edge is the lowest edge of one of its 
extremities in any graph invariant by 7e- On the contrary, J,"^=J, 
transforms an arbitrary graph into a graph invariant by 7e. 

The pruning J,^ keeps for each node / the adjacent edges which are 
followed by a second couple of minimal weight. These edges are those 
linking /' with one of its lowest neighboring nodes. 
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A morphological characterization of k-steep graphs 
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Eroding k-steep graphs 



Consider a k-steep graph J,'' G = G 
erosion eG = (£ee,enn), where £e = 
e{i)G = eG and e('")G = ee^^'-^'G. 



= (e, n), with /c > 2. We define the 

^en^ne Sid — ^ne^en- 



It does not change the catchment basins if we assign to the regional 
minima the weight 0, whereas all other nodes and edges have weights > 
(like that as soon the erosion assigns the value to a node or an edge, 
they remain equal to for all subsequent erosions). 



cr 
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Eroding k-steep graphs 



Consider a k-steep lexicographic track t of a flooding graph G made of a 
series of non increasing flooding pairs (/'i, ei), (/2, 62), ... , {i^, e/^). Node 
and edge of each flooding pair have the same weights. As t has a minimal 
lexicographic weight, each of its flooding pairs, except the first one 
constitutes one of the lowest adjacent flooding pair of the previous flooding 
pair. For this reason the erosion eG assigns to the edge e/, the weight of 
the adjacent edge e/,+i and to the node /'/, the weight of the adjacent node 
ih+i- In other workds successive erosions (Cee, £„/!) let glide the value of 
each pair upwards in the track. If this track is of length k, after k — 1 
erosion, the weight of the ultimate pair will have reached the first one. 
If such a track is of length / < /c, it ends with a pair in a regional 
minimum, with the value 0. Successive erosions (£ee,£n/i) also let glide the 
value of each pair upwards in the track and the last pair, with value also 
moves upwards and reaches the pair (/', y) after / — 1 erosions. During the 
next erosions, the value of the pair (/'i, ei) remains stable and equal to 0. 
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Eroding k-steep graphs 



Consider a k-steep flooding graph G =1^^ G. During the k — 1 successive 
erosions, the values of the flooding pairs glide upwards along the k-steep 
lexicographic path. For the erosion / < k, the (/ -|- l)th pair (s, st) has 
reached the pair (/', //'). Hence the edge ij remains one of the lowest 
adjacent edge of the node / and both share identical weights. 
But (s, St) belongs to the flooding graph and verifies = £nes{s) and 
/Is = enes{s) and so does (/, ij). 

This shows that the graph e^'^G still is a flooding graph. 



Theorem 



For an ordinary flooding grapfi G , G is a k-steep flooding grapli and for 
I < k, eC) I'' G still is a flooding graph. 
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Eroding flooding graphs 



Consider a flooding graph G = (e, n). The erosion £„ = Cne^en enlarges 
the regional minima. Hence if G is invariant by cpn it is still the case for 
the graph (e, 

On the contrary, after the erosion £e = Sen^ne, the graph {See, n) is not 
necessarily invariant for 7e. One has to prune the edges which are not the 
lowest edges of one of their extremities: \. egfi. 

Repeating the operator ^G = (J, £ee,£nn) produces a decreasing series of 
partial graphs ^(")G = ^'-"-^^ G we will now characterize. 
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Two equivalent modes of pruning. 



Theorem 



For m > 1, and defining = identity, ttie operators ^ applied to 

G and |'"+^ applied to I"" G keep and discard the same edges. The 
operators ^('^^ and l^^^ produce partial spanning graphs with the same 
edges and nodes. 

Proof: We prove it by induction. 

a) m = 1 : ^ and J,^ select the edges linking a node with its lowest node 
(this covers the case where the edge belongs to a regional minimum) 
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Two equivalent modes of pruning. 



b) We suppose that ^('"^^)g and J,'" G have the same edges and show 
that it is still the case for m = m + 1. 

Consider a couple (/', ij) of ^("'^^)g. It belongs to a m — 1-steepest track. 
The second pair (JJk) of this track also belongs to a [m — 1) —steepest 
track and holds the weight of the lowest pair of the track. The operator ^ 
applied to ^('"^i)G assigns this weight to the edge ij but not necessarily 
to the node / if there exists another pair (/, /'/) with a lower weight after 
this erosion. In this case the edge ij is discarded by the operator ^. But it 
is also discarded by the operator l'"^^ G, as no (m + l)-steepest path 
passes through ij. 
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Constructing k-steepest graphs 



The operator J,'" G is of theoretical interest but of poor practical value, as 
it is based on a neighborhood of size m. On the contrary the operator ^ is 
purely local and uses a neighborhood of size 1 : an erosion from node to 
node, an erosion from edge to edge and the suppression of any edge which 
is not the smallest edge of one of its extremities. 

Repeating m times the operator ^ produces a graph which has the same 
edges as the operator G. It is then sufficient to restore the original 

weights of edges and nodes of the graph G onto the graph ^("'^G to 
obtain the same result as G. 



r 1 II I mill I II I II II 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigl 2012 February 29 107 / 201 



Characterizing k-steepest graphs 



Theorem 



A graph G is a k-steepest graph if and only if for each I < k , e^'^G is a \ 
flooding graph J 

Proof: If a graph G is a k-steepest graph, then G =\.^ G and we have 
already established that for / < k, e''^ J,'' G is a flooding graph. 
Inversely suppose that for each I < k , e^'^G is a flooding graph, then 
e^'^G has the same edges as G. As a matter of fact e^'^G is invariant by 7e 
and each edge is the lowest edge of a node: J, e^'^G = e^'^G. 
But eC) G = ££('-1) G and ; ee^'-'^^ G = {I eee, enn)e^'-'^^ G = ^£('-^) G 
(considering G = e^^^G). It follows that for / < /(, we have £(')G = ^^G. 
The operator ^ applied to /c times to G never suppresses an edge from G. 
But we know that ^^''^-^'G and G have the same nodes and edges, 
showing that indeed G is a k-steepest graph. 
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Deriving an algorithm for delineating the catchment basins 
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The catchment basins of k-steepest paths 



It is now possible to imagine an algorithm associated to the operator ^. If 
G is a flooding graph, we assign a a distinct label to each regional 
minimum and a weight 0. The operator ^ propagates the weights upwards 
along the steepest tracks as it is repeated. In parallel, we also propagate 
the labels: every time a pair (/', ij) takes the weight from a pair 
containing a labeled node j, the node / takes the label of j. Like that, as 
the labeled zones with weight expand, their labels also expand. 
We now have to consider two cases. In the first case, we construct 
restricted catchment basins separated by watershed zones. In the second, 
we create a partition by assigning each node to one and only one basin, at 
the price of arbitrary choices. 
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Catchment basins and watershed zones 



During the successive operations ^, every time a pair (;', ij) takes the 
weight from a pair containing the node j, the node / takes the label of J 
provided the pair (/', ij) is unique. If there are two equivalent pairs (;', ij) 
and (/', //), it means that there exist two minimal lexicographic tracks 
starting at / towards one or two distinct minima. If j and k hold the same 
label, this label is assigned to /'. If on the contrary they are distinct, we 
have 2 possibilities: 

• we assign to / a label Z, indicating that it belongs to a watershed 
zone. The upstream of / also will get this same label Z. The regions 
with label Z belong to the watershed zone and the other labeled 
regions are restricted catchment basins. 

• we assign to / one of the labels (either randomly, or by applying an 
additional rule for braking the ties), producing a partition in 
catchment basins with an empty watershed zone. The result is not 
unique and depends on the succession of choices which have been 
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Creating a partition 




On the left a flooding graph where the minima have weights equal to 
and their labels are indicated by distinct colors. The next two figures show 
the propagation of the weights and of the labels as the operator ^ has 
been applied twice in a sequence. 
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Partial conclusion 



Starting with a node or edge weighted graph we have shown how to 
extract from it a flooding graph where nodes and edges are weighted. 
The operators J,"' G extract from the graph G partial graphs which are 
steeper and steeper flooding graphs. The series of partial graphs J,™ G is 
decreasing with m. 

The operator ^ permits an iterative construction of J,'" G, using only small 
local neighborhood transformations. 
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The scissor operator and the watershed partitions 

In the previous section we introduced pruning operators which extract 
from a flooding graph k-steep flooding graphs. These operators do not 
make any choice among the k-steep flooding graphs, they take them all. 

We now introduce operators aiming at creating partitions : they extract 
minimum spanning forests from the flooding graph by pruning. Contrarily 
to the preceding operators, they do make arbitrary choices among 
equivalent edges to be suppressed. 
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The flooding pairs 



The previous section has shown how to extract from a node or edge 
weighted graph partial graphs which are flooding graphs. In a flooding 
graph, each node has at least one adjacent edge with the same weight. We 
consider here the nodes and edges outside the regional minima. There 
exists at least one (in general several) one to one correspondance between 
each node outside a regional minimum and one of its adjacent edges with 
the same weight. For each such one to one correspondance, we call 
flooding pair, the couple of node and edge which have been associated. 
They hold the same weight. 
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The flooding pairs 



Let us show how to construct such a one to one correspondance. If an 
edge is the lowest adjacent edge of only one of its extremities, then they 
form a flooding pair ; this is in particular the case if the other extremity of 
the edge has a lower weight. 

If on the contrary, it has the same weight as both its extremities, it 
belongs to a plateau. As this plateau is not a regional minimum, there 
exists a pair of neighboring nodes, a node s inside the plateau and a lower 
node t, outside. This edge (s, t) forms a flooding pair with the node s. As 
the plateau is connected, there exists a tree spanning its nodes, having s 
as root. There exists a unique path between s and each node / of the 
plateau. The last edge on this path before reaching / and /' itself form a 
flooding pair. 
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The scissor operator creates a partition of the nodes 



We have shown that there exists at least one (in general several) one to 
one correspondance between each node outside a regional minimum and 
one of its adjacent edges with the same weight. For each such one to one 
correspondance, we call flooding pair, the couple of node and edge which 
have been associated. 

To each such one to one correpondance we associate a a pruning operator 
X, called scissor. This operator suppresses all edges outside a regional 
minimum which do not form a flooding pair with one of their extremities. 

After applying the operator ;t to a flooding graph, there exists one and 
only one path from each node to a regional minimum: the catchment 
basins of the minima partition the nodes, and the watershed zones are 
empty 
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The drainage minimum spanning forest 
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The drainage minimum spanning forest 



After contracting all edges in the regional minima of a flooding graph and 
applying the scissor operator x g^ts a spanning forest, where each tree 
is rooted in a minimum: 

• the resulting graph is a partition where each connected component 
contains a regional minimum node : there exists a path linking each 
node outside the minima with a minimum. 

a each connected component is a tree as the number of nodes is equal 
to the number of edges (the number of flooding pairs) plus one (the 
regional minimum node). 

a it is a minimum spanning forest : each node is linked to its tree 
through one of its lowest neighboring edge. The total weight of each 
tree is thus equal to the total sum of the nodes outside the regional 
minima. This total weight is independent of the particular choice 
made by 

' r 1 ill I mill I II I II II I 



Fernand MeyerCentre de Morphologie Mathe IVaters/iec/s, waterfalls, on edge or node weigl 2012 February 29 119 / 201 



The drainage minimum spanning forest 



Expanding again the regional minima and replacing them by a MST of 
each regional minimum creates again a minimum spanning forest, identical 
to the preceding one outside the minima. 
Its total weight is computed as follows: 

• if Mi is a regional minimum with n nodes of weight A;, the weight of 
its MST is equal to {n — 1) * A/ 

• each other node contributes by its weight, which is also the weight of 
its lowest adjacent edge. 

Each particular forest is based on a particular scissor operator x and of a 
particular MST in each regional minimum. We call (p the operator which 
extracts from a flooding graph such a MSF (minimum spanning forest). 
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The catchment basins 



G = a flooding graph, C = a minimum spanning forest: 

• Each node belongs to a regional minimum or is the origin of a unique 
never ascending path leading to a regional minimum 

• each tree of C contains a unique regional minimum ; its nodes form 
the catchment basin of this minimum. 
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Catchment basins of increasing steepness 



G = a flooding graph, J,'" G still is a flooding graph and Cm = (p ■l'" G a 
minimum spanning forest of steepness m. 

As for m > / : J,"" G Cj,' G , any minimum spanning forest Cm = (p ^."^ G 
is also a minimum spanning forest of J,' G. 

For increasing values of m, the number of forests of steepness m decreases. 
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Illustration 
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A flooding graph associated to a distance function 



We chose as topographic surface the distance function expressed on the 
nodes of a hexagonal grid to two binary connected sets, encoded with the 
value 0. The edge weights are obtained by the dilation Sen- Like that the 
lowest neighboring edges of a node connects it with its neighbors with 
equal or lower weights. In the following figure, the edges with the same 
weight have the same color (cyan = 4, magenta = 3, green = 2). 

3 — 2 — 2 
/ \ / \ / \ 
4—3—3—3 
\ / \ / \ / 

4—4—4 
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The flooding graph of a distance function 



As the minima have 2 nodes each, the pixel graph is invariant by (pn- The 
dilation Sen assigns weights to the edges and creates a flooding graph. 
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The minimum spanning forest by keeping one lowest 
neighboring edge for each node 



The scissor x leaves one lowest neighboring edge for each node. There 
exists a huge number of choices for X- The following illustration shows a 
particular scissor x producing an unexpected partition in two catchment 
basins (the catchment basins should be the half plane separated by the 
mediatrix of both binary sets, as illustrated in the previous slide). 
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Two geodesies for the flooding ultrametric distance, leading to an 
unexpected partition. 
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Looking one node further 



The operator J,^ G leaves only the edges linking a node to its lowest 
neighbors. The following figure shows the remaining edges. The green zone 
represents a restricted catchment basin. The yellow zone is an extended 
catchment basin containing the watershed zone. 
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The minimum spanning forest by keeping one edge towards 
a lowest neighboring node 



The pruning x. J|= X ^ leaves for each node one edge towards a lowest 
neighboring node. The following figure shows one such solution and the 
resulting partition. 
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The catchment basins, skeletons by zone of influence for 
lexicographic distance functions 
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The geodesies of the k-steepest graphs 



The operator J,'* G prunes the flooding graph and leaves only NAPs with a 
steepness equal to k. From each node of the graph starts a NAP whose k 
first edges have a minimal lexicographic weight. If one follows such a path, 
the next k edges following each node as one goes downwards along the 
path also has a minimal lexicographic weight. This shows that each NAP 
is a geodesic line for a lexicographic distance function lexdist/< we define 
below. 

After the pruning J,'^ some nodes belong to two or more catchment basins. 
In order to obtain a partition, one applies the scissor operator x leaving for 
each node only one lowest adjacent edge. Like that the thick watershed 
zones are suppressed and the catchment basins form a partition. 
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Constructing a watershed partition as the skeleton by 
zones of influence of the minima 



It is possible to obtain the same result, if one labels the regional minima 
and computes for the other nodes the shortest lexicographic distance 
lexdist/t to the minima. If one applies a greedy algorithm, one may in 
addition propagate the labels of the minima all along the geodesies and 
construct a partition of the space. If a node is at the same lexicographic 
distance of two nodes, a greedy algorithm will arbitrarily assign to it one of 
the labels of the minima. This is quite similar to the operator x which also 
does arbitrary choices. 
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Defining a lexicographic distance along non ascending 
paths. 



Consider a flooding graph. In a NAP, each node except the last one forms 
with the following edge a flooding pair, i.e. they have the same weights. 
And these weights are decreasing as one follows the NAP downwards. The 
k first values, starting from the top, may be considered as a lexicographic 
distance of depth k. In what follows we give a precise definition of such 
distances. 

The shortest distances and their geodesies may be computed wih classical 
algorithms. Propagating the labels of the minima along the geodesies 
during their construction constructs the zones of influence of the minima, 
i.e a partition into catchment basins. The solution is not necessarily 
unique. However the number of solutions decreases with the depth of the 
lexicographic distance which is considered. 
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Lexicographic distances of depth k 
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Comparing NAPs with the lexicographic order. 



Let S be the set of sequences of edge or node weights, i.e. elements of T . 
Let Sk be the set of sequences with a maximal number k of elements. For 
a sequence s G Sk, we define the lexicographic weight Wk{s) : it is equal 
to 00 if s is not a NAP and equal to s itself otherwise. 
We define an operator first/< which keeps the k first edges and nodes of 
any NAP, or the NAP completely if its length is smaller than k. The 
operator ikstk maps any sequence of S into Sk- 

We define on the NAPs of Sk the usual lexicographic order relation, which 

we will note ^, such that: (Ai, A2 \k) ~< 1^2 I^k) if either 

Ai < ^1 or A/ = }ii until rank s and Aj+i < /^s+i- We define a < b as 
a ~< b or a = b. 
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Comparing NAPs with the lexicographic order. 



Like that, it is possible to compare any two sequences si and S2 of Sk by 
comparing their weights: 

• if si and S2 are not NAPs, then Wk[si) = m//<(s2) = oo and we 
consider that si = S2 (they are equivalent) 

• if one of them, say si, is a NAP and not the other, then 

Wk{si) -< Wk{s2) = 00 and si -< S2 

• if both of them are NAP, then they compare as their weights: 

si S2<^ {wk{si) -< Wk{s2)} and si = S2 <^ {m//<(si) = Wfc(s2)} 

If si and 52 are NAPs belonging to S, we compare them with: si ~< S2 
Wk [first/,(si)] ^ Wk [firstfc(s2)] 
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Defining an "addition" operator, comparing sequences 



For NAPs of S we define the operator ffl/( called "addition", which 

operates as a minimum: 

m , a if a ^ 6 , , , r- 
3Skb= , , Z Va, fa e S 

D if D ^ a 

The ffl/f operation is associative, commutative and has a neutral element 
00 called the zero element: a E/^ oo = oo S/^ a = a 



cr 
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Defining a "multiply" operator, concatenating sequences 



The operator M^, called "multiplication", permits to compute the 
lexicographic length A[A) of a sequence, obtained by the concatenation of 

two sequences. Let a = (Ai, A2 /^k) and b = ^2 ^k) we will 

define a^k b by: 

• if a or ifa is not a NAP then aMk b = 00 

• if a and b are NAPs and Ak < }ii then aMk b = ca 

• if a and b are NAPs and \k > fii then a Kl/c i> = [first/((a O b)] , 
where a > 6 is the concatenation of both sequences. 

In particular a 00 = 00 a = 00. so that the zero element is an 
absorbing element for lEj/^. 
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Algebraic shortest paths algorithms 
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A path algebra on a dioTd structure 



The operator Mj^ is associative and has a neutral element called unit 
element: a^i^O = a. The multiplication is distributive with respect to the 
addition both to the left and to the right. 

The structure (S, ffl/<,KI/() forms a dio'i'd, on which Gondran and Minoux 
defined a path algebra [TU], where shortest paths algorithms are expressed 
as solutions of linear systems. 
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Transposing the dioTd to square matrices 



Addition and multiplication of square matrices of size n derives from the 
laws ffl and : for A = (a/,-), B = {bij), ij E [1, n] : 

C = AmB = (cij) O Cij = aij ffl bij V/, j 
C = AmkB = (dj) ^Cij= E aik bkj ^ij 

l<k<n 



where ^ is the sum relative to ffl. 

ffl 

As there is no ambiguity, for ^ a//( 6/^-, we simply write ^ ajkb^j 

l<k<n l<k<n 
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Unity and zero matrices 



With these two laws, the square matrices also become a dioVd with 

£ £ 



zero matrix £ 



and unity matrix E 



We write = E 



e £ 

..e 

....e 

e 

e.... 

£ e 



I jl Ij I II II II I II n 



II I llll I 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigl 2012 February 29 142 / 201 



Lexicographic length of paths in a graph 



If G = [X, U] is a weighted graph with 

• a set X of nodes, numbered / = 1, , N. 

9 a set U of edges u = (/,_/') with weights s,y G S. 
The incidence matrix A = (ay) of the graph is given by: 



SI 



(ij) e u 

sinon 
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A path algebra on a dioTd structure 



To each path = (/'i, 72, ....//<) of the graph, one associated its 
k-lexicographic weight w{}i) = s;^/^ Kl/^ s/j/j ■■■■^k ^ik-vk' which is 
different from oo only if the path is a NAP. 

If 71 is a never increasing sequence, then Wf^ [first/<(7r)] = first/((7r) 



Shortest paths in the graph 



The shortest paths for the lexicographic distance between any couple of 
nodes may be com puted thanks to A" or /\(") = E ffl /\i ffl g 



Lemma 



where: 

• C," is the family of paths between / and _/, containing n + 1 nodes 
(not necessarily distinct) 

• c\"^ is the family of paths between /' and _/', containing at most n + 1 
nodes (not necessarily distinct) 

Defining /\' = E E /\, as ffl est idempotent (a E a = a Va G S), we 
have: = A^"^ 
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A path algebra on a dioTd structure 



In a graph with N nodes, an elementary path has at most N nodes, 
separated by A/ - 1 edges. Hence, necessarily Z^^'^) = /\('^"^) and A*, the 
limit of /\(") for increasing n, is also equal to A^'^ ^\ 
Thanks to A*, the computation of the catchment basins is immediate. If 

mi is a regional minimum, then a node / belongs to its catchment basin if 

EB 

andonlyif/\;^^< E AJ^^ 

A* may be obtained by the successive multiplications : 

A, A'^ = AmkA, /^4 = ^2^^^2^ ^2- = 42--! ^2-1 ^ntil 

2' > N-1,, i.e. / > log (A/ - 1) 
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A path algebra on a dioTd structure 



A* verifies A* = EmAmkA*. 

Multiplying by a matrix B: A* B = BmAMkA*B. Defining Y = A* B 

shows that A* B is solution of the equation 

Y = BS A^k Y ■ Furthermore, it is the smallest solution. 

With varying B it is possible to solve all types of shortest distances: 

• The smallest solution Y = EmAMkY yields A* E = A* 



» with 6 



one gets A* B, the i-th column of the matrix A* 



that is the distance of all nodes to the node /'. 
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Solving linear systems 



Gondran and Minoux have shown that most of the classical algorithms 
solving systems of linear equations (Gauss, Gauss-Seidel, etc.) are still 
valid in this context and correspond often to known shortest paths 
algorithms defined on graphs. We now give a few examples. 



II II II I II II 



I 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigl 



12 February 29 148 / 201 



The Jacobi algorithm 



Setting B 



and 



the iteration 



YW =A* yC'-i) ffl B converges to Y^^^ = A* * B 
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The Gauss Seidel algorithm 



A is decomposed as 4 = /. ffl E ffl (7, where L is an inferior triangular 
matrix, E the unity matrix E, and U an upper triangular matrix. The 
upper part of L and lower part of U have the value e. 

The solution ofy = /4*yE6is obtained by the iteration : 

yik) ^ iY{k-\) g ^Y(k) g g 

This algorithm is faster as Jacobi's algorithm, as the product UY^'^^ 
already uses intermediate results, freshly computed during the current 
iteration Y'^''\ 
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The Jordan algorithm 



The Jordan algorithm is used in classical linear algebra for inverting 
matrices, by successive pivoting. In our case, where the shortest paths are 
elementary path the algorithms is: 
For k from 1 to A/ : 

For each / and J from / to N, do: aij = ajj ffl aik * a^j 
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The greedy algorithm of Moore-Dijkstra 



Gondran established the algebraic counterpart of the famous shortest path 
algorithm by Moore-Dijkstra. 

Theorem (Gondran): Let Y = A* B be the solution Y = AYS B, for 

ffl 

an arbitrary matrix B. There exists then an index /'o such that = J^b;. 
The smallest b is such solution : = b/^ 

Each element of V = AY ffl 6 is computed by 

ffl ffl 
y/c = E 3kj *yjSbk = £ a/y- * yj E aki^yi^ E bk 

Suppressing the line and column of rank /q and taking for b the vector 

b^k^ = akioYio S bk, one gets a new system of size A/ — 1 to solve. 

The Moore Dijkstra algorithm can also directly be computed on a flooding 

graph G, as presented below. 
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Shortest paths algorithms on the graph 



• The shortest path algorithm of Dijkstra is first presented [22]- We 
show that for a lexicographic distance of depth 1, it becomes 
algorithm for constructing the minimum spanning forest of Prim. 

• The core expanding algorithms take advantage of the particular 
structure of the lexicographic distances. They are faster than the 
Dijkstra algorithm and better suited to hardware implementations. 
For a lexicographic distance of depth 2, they produce the same 
geodesies as the topographic distance. 
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Architecture of shortest path algorithms. 



The shortest path algorithms below are applied on a graph which is 
invariant by the opening 7e (the regional minima may ad libitum be 
contracted beforehand). A domain D is used and expanded, containing at 
each stage of the algorithms the nodes for which the shortest distance to 
the minima is known. Initially the minima are labeled and put in D with a 
value 0. We say that a flooding pair (J,JI) is on the boundary of the 
domain D, if j is outside D and / inside D. In this case we say that "j 
floods /" , or "/ is flooded by /' . We say that j belongs to the outside 
boundary 9+D of D and / to the inside boundary 9^D of D. 
The domain D is progressively expanded by progressive incorporation of 
nodes in 3+D belonging to flooding pairs on the boundary of D. 
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The Moore Dijkstra algorithm 



The shortest path algorithm by Moore-Dijkstra constructs the distances of 
all nodes to the minima in a greedy manner. It uses a domain D which 
contains at each stage of the algorithm the nodes /' for which the shortest 
distance S*^{i) and label A(/) is known. 
Initialisation: 

The nodes of the regional minima are labeled and put in the domain D. 
Repeat until the domain D contains all nodes: 

For each flooding pair (J, (Jl) on the boundary of D, estimate the shortest 
path as S,{j) = ej,^kSl{l). 

The node with the lowest estimate is correctly estimated. If the 
corresponding flooding pair is (s, st) : 

• D = DU{s} 

• Sl{s) = S,{s) 

• A(s)=A(t) 
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Correctness of the Moore Dijkstra algorithm 



The node with the lowest estimate is correctly estimated and is introduced 
in the domain D. If is necessarily the shortest path, as any other path 
would have to leave D through another boundary flooding pair with a 
higher estimate. 



Controlling the Moore Dijkstra algorithm 



The Moore Dijkstra may be advantageously controlled by a hierarchical 
queue structure. Each node, as it gets is estimate is put into the HQ. As 
long the HQ is not empty, the extracted node is among the nodes with the 
lowest estimate, the one which has been introduced in the HQ first. 
The HQ has thus the advantage to correctly sequence the treatment in the 
presence of plateaus: it treats the nodes in the plateaus from their lower 
boundary inwards. The processing order is proportional to the distance of 
each node to the lower boundary of the plateau. 
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The Moore Dijkstra algorithm : case where k = 1 



The algorithms remains the same but the computations are simplified as 
Sk(J) = eji Kli is simply the weight of the node and the edge in the 
flooding pair {J, (Jl). In other terms, the domain D is expanded by 
introducing into D the flooding pair on the boundary of D with the lowest 
weight. This corresponds exactly to the algorithm of PRIM for 
constructing minimum spanning forests. The same algorithm has been 
used in [IB] for constructing the waterfall hierarchy. 
This is not surprising, as for k = 1, the lexicographic weight of a NAP 
simply is the weight of the first edge. The distance is in this case the 
ultrametric flooding distance. 

Remark: There is a complete freedom in the choice of the flooding pairs 
which are introduced at any time into D. A huge number of solutions are 
compatible with this distance, some of them quite unexpected as 
illustrated by an example given above. 
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The core expanding shortest distance algorithm 



Due to the particular structure of lexicographic distances, it is possible to 
identify another type of nodes for which the shortest distance may 
immediately be computed. Let be the set of nodes in the boundary 

D with the lowest valuation first/(_i(cyj^). For each t G 9 and for 

each s G flooding t we do Cfs fii'st^-i ^l{t) producing a NAP of 
length k. The value of Sl{t) inst^-i is simply obtained by appending the 
weight of s to first/(_i ^l{t). As t belongs to 9 0, this value is the 
smallest possible. 

This analysis shows that each node of permits to introduce into D, all 
its neighbors by which it is flooded, whatever their weight. This algorithm 
is more "active" as Dijkstra's algorithm, as each node inside is able to 
label and introduce into D all its flooding neighbors at the same time. 
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The core expanding shortest distance algorithm 



The algorithm is the following; 
Initialisation: 

The nodes of the regional minima are labeled and put in the domain D. 
Repeat until the domain D contains all nodes: 

Let D be the subdomain of D of nodes flooded by nodes outside D 



with the lowest valuation first/(_i(c/^). For each t G ^'^'^ ^ flooded 
by G a+D : 

• D = DU{s} 

• A(s)=A(t) 




Controlling the core expanding shortest distance algorithm 



The core expanding shortest distance algorithm may also be 
advantageously controlled by a hierarchical queue structure. Each node, as 
it is introduced in D is put into the HQ. As long the HQ is not empty, the 
extracted node is among the nodes with the lowest estimate, the one 
which has been introduced in the HQ first. This node gives its label and 
the correct distances to all its neighbors outside D by which it is flooded. 
If a node extracted from the HQ belongs to a plateau, as it has been 
introduced in the HQ, it is closer to the lower boundary of the plateau 
than other nodes which may have been introduced in the HQ later. 
The algorithm is particularly suitable for a hardware implementation: each 
node is entered in the HQ only once. On the contrary, with the 
Moore-Dijkstra algorithm a node may be introduced several times in the 
HQ, as its estimate may vary before it gets its definitive value. 
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The core expanding shortest distance algorithm: case 
where k = 1 



The domain has been defined as the subdomain of D of nodes 

flooded by nodes outside D with the lowest valuation first/c_i(c/^). For 
/c = 1, we have first/<_i(c/^) is empty, and d D occupies the whole domain 

d D. This means that any node in 9^D can be expanded by appending 
one of the outside node through which it is flooded. The algorithms for 
constructing graph cuts by J.Cousty find also their place in this context [2] 
We have illustrated this situation earlier, showing that the minimum 
spanning forests with steepness may be absolutely unexpected. 
Controling the algorithm with a hierarchical queue limits the anarchy to 
some extend. 



r 1 ill mill I II I II II I 



Fernand MeyerCentre de Morphologie Mathe Watersheds, waterfalls, on edge or node weigl 2012 February 29 162 / 201 



The core expanding shortest distance algorithm: case 
where k = 2 



\f k = 2, then the valuation first/<_i(c/^'/)) = firsti(c/|(/)) is the valuation 
of the node / itself. This means that contains the nodes with lowest 

weight belonging to 9^D. If / is such a node, it introduces into D each of 
its neighbors j belonging to 9+D, each with a valuation firsti(cy|(j)) equal 
to its weight nj. 

If in addition, one uses a HQ for controlling the process, we get the 
classical algorithm for constructing catchment basins |17].|19|. 
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The classical watershed algorithm 



Label the nodes of the minima and them in the domain D, each with a 
priority with a weight. 

As long as the HQ is not empty, extract the node j with the highest 
priority from the HQ: 

For each unlabeled neighboring (on the flooding graph) node / of J : 

* label{i) = label{j) 

* put / in the queue with priority i// 

As a matter of fact, this algorithm has first been derived from the 
watershed line, as zone of influence of the minima for the topographic 
distance, defined below. 
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The topographic distance 



Consider an arbitrary path n = (xi,X2, ...,Xp) between two nodes xi and 
Xn. The weight Vp at node Xp can be written: 

Vp = Vp- l/p_i + l/p_i - l/p_2 + Vp-2 - Vp-3 + ....+ 1/2 - t^l + 

The node /c — 1 is not necessarily the lowest node of node k, therefore 

Vk-l > £nVk and Vk - Vk-l <Vk- tnVk- 

Replacing each increment Vi^ — by Vi^ — En^k will produce a sum 
Vp — £nVp + i/p_i — £nVp-i + .... + 1/2 — enV2 + Vi which is larger than 
Vp. It is called the topographic length of the path n= (xi, X2, Xp). The 
path with the shortest topographic length between two nodes is called the 
topographic distance between these nodes. 
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The topographic distance 



The path with the shortest topographic length between two nodes is called 
the topographic distance between these nodes [23], [12]. It will only be 
equal to Vp in the case where the path (xi,X2, ...,Xp) precisely is a path of 
steepest descent, from each node to its lowest neighbor. In other terms the 
topographic distance and the distance Iexdist2 have the same geodesies. 
If we define the toll to pay along a path n = (xi,X2, ...,Xp) as Vi for the 
first node and V; — £„!// for the others, then the lowest toll distance for a 
node Xp to a regional minimum will be Vp if there exists a path of steepest 
descent from Xp to xi. In other terms, Xp and xi belong to the same 
catchment basins if one considers the topographic distance. But they also 
belong to the same catchment basin if one considers the depth 2 
lexicographic distance, as, by construction, x/^ is the lowest neighbor of 

Xk-l- 
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Distance on a node weighted graph based on the cost for 
travelling along the cheapest path 



We recall the toll distance presented earlier. 

The weights are assigned to the nodes and not to the edges. Each node 

may be considered as a town where a toll has to be paid. 

Cost of a path: The cost of a path is equal to the sum of the tolls to be 

paid in all towns encountered along the path (including or not one or both 

ends). 

Cost between two nodes: The cost for reaching node y from node x is 
equal to the minimal cost of all paths between x and y. We write 
tolldist(x, y). If there is no path between them, the cost is equal to oo. 
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I 



ustration of the cheapest path 































^X2) 










— <3) 



In this figure, the shortest chain between x and y is in red and the total 
toll to pay is 1 + 1 + 2 + 2 = 6 
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Reconstruction of an image by integration 



Finally, any image f may be considered as the global toll of its pixel graph 
if one takes: 

• as reference nodes the regional minima of the image. Each of them 
has as toll its altitude. 

« as local toll for all other nodes, the difference between their altitude 
and the altitude or their lowest neighbor: g = f — ef 

If in addition, we give a different label to each regional minimum, we may 
as previously propagate this label along each smallest toll path. We obtain 
like that a tessellation: to each minimum is ascribed a catchment basin: 
the set of all nodes which are closer to this minimum then to any other 
minimum. 
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Inversely: the catchment basins of the cheapest paths for a 
distribution of tolls. 



Inversely we may chose a number of starting nodes called roots, with a toll 
to pay and for all other nodes, the toll to pay for reaching or crossing this 
node. The toll to pay constitutes a topographic surface where each root is 
a regional minimum. In addition we propagate the labels of the minima 
along the geodesies of the cheapest distance, and get a partition of the 
nodes. As a result one get a partition of the space, where each region with 
a given label is the catchment basin of the root with the same label for 
this distance function. 
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Assigning to each node the global toll for reaching it, 
starting from one of the roots. 
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ustration of the cheapest path 



For each node we have indicated the local toll value (left value) and the 
global toll to pay to reach the closest reference node, who shares the same 
color (label): 

• Each reference node became a regional minimum of the graph. 

• The local toll of any other node is equal to the difference between its 
global toll and the global toll of its lowest neighbor. 

• Each non reference node got its value and label from one of its lowest 
neighbors. 

• The values of the nodes are computed and the labels propagated 
along a path of steepest descent. 
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I 



ustration of the cheapest path 




The catchment basins of this surface are the SKIZ of the minima, both for 
the topographic distance and for the depth 2 lexicographic distance. Each 
node is the extremity of a geodesic line, which is the same both for the toll 
distance and for the lexicographic distance of depth 2 computed on the 
same topographic surface. 
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Top down or bottom up 
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The influence of the depth k 



For increasing values of k, the domain 9 D becomes smaller and smaller, 
indicating that the number of equivalent catchment basins compatible 
with a given lexicographic depth is reduced as the value of k becomes 
bigger. This is in accordance with the fact that with increasing values of 
k, the pruning J,'' becomes more and more severe. 
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Obtaining catchment basins with k — steepness 



After applying the operator J,'' on a flooding graph G, there remain only 
NAP with k steepness. 

The same is true if we apply the operator ^^'^^^^ in order to prune edges of 
G and subsequently restore the initial weights of the edges onto the 
remaining edges. 

If the only NAPs remaining in the graph have a k steepness, any shortest 
path algorithm with a lower steepness will extract them In particular the 
most simple algorithms for the distances di or c/2 will extract catchment 
basins of steepness k. 
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ustration : lexicographic distances 



The follownig 3 images show respectively the shortest lexicographic 
distances of depth 1, 2 and 3. If a path is the shortest path for a 
lexicographic distance of depth k, it also is a shortest path for a 
lexicographic distance of smaller depth. The following three figures present 
the lexicographic distances of depth 1, 2 and 3. The partition of 
catchment basins for the distance 3 is also solution for the distance 2 and 
1. Similarly the partition for distance 2 is also solution for distance 1. The 
number of solutions decreases with the lexicographic depth. 
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I 



ustration : lexicographic distance of depth 1 
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I 



ustration : lexicographic distance of depth 2 




I 



ustration : lexicographic distance of depth 3 




Erosion and pruning 



Repeating the operator = (J, eee,enn) produces a decreasing series of 
partial graphs ^("^G = ^^("^-"-'G, which are steeper and steeper. In the 
following figures, we present in red the edge which is not the lowest edge 
of one of its extremities. After pruning this edge, the operator is applied 
again. Applying ^ a number n of times is equivalent to constructing the 
partitions compatible with a SKIZ for a lexicographic distance of depth 
n + 1. In our case, ^'^^G produces a graph with the same edges as the 
pruning x applied to the graph where each edge has been weighted by its 
lexicographic distance of depth 4 to the nearest minimum, illustrated in 
the previous figure. 
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Erosion and pruning 
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The hierarchy of nested catchment basins 
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Watershed and waterfalls 



The waterfall hierarchy has been introduced by S.Beucher in order to 
obtain a multiscale segmentation of an image [2],E],E]. Given a 
topographic surface Si, typically a gradient image of the image to 
segment, a first watershed transform produces a first partition tti. 
The waterfall flooding of Si consists in flooding each catchment basin up 
to its lowest neighboring pass point, producing like that a topographic 
surface with less minima. The watershed segmentation of this surface 
produces a coarser partition 1x2 where each region is the union of a 
number of regions of TZi. 
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Chaining the waterfall floodings 



Flooding each catchment basin of Si up to its lowest pass point produces 
a new topographic surface S2 which will be submitted to the same 
treatment as the initial surface Si. Its watershed transform produces a 
second partition 712- The partition 712 is coarser than tti as each of its tile 
is a union of tiles of Tii.The following figure shows how the flooding of Si 
produces S2, which, in the second figure, has also been flooded. 



The same process can be repeated several times until a completely flat 
surface is created. The partitions obtained by the watershed construction 
on the successive waterfall floodings are coarser and coarser : they form a 



hierarchy. 
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The waterfall hierarchy 



A 1 dimensional topographic surface is represented through the altitude of 
its pass points in the figure below. The first flooding assigns weights to the 
nodes and its watershed construction produces 4 catchment basins, 
separated by 3 watershed lines. The second flooding has only 2 catchment 
basins separated by 1 watershed line. The last image orders the watershed 
lines of the initial image into 3 categories, in cyan the watershed lines 
which disappeared during the first flooding, in dark blue those which 
disappeared after the second flooding and in red the one which survived 
the second flooding. 
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Watershed and waterfalls 



Let us come back to the watershed on weighted graphs. The watershed of 
the topographic surface produces a partition Tii into catchment basins, 
represented by its region adjacency graph RAGi. The first flooding floods 
each catchment basin up to its lowest neighboring pass point. This 
corresponds to the erosion £ne of the graph RAGi. The theory of the 
watershed on weighted graphs can now be applied on this graph. The 
resulting watershed appears in form of minimum spanning forest MSFw 
each tree of the forest spans a catchment basin of the partition tc2- 
The next level of the hierarchy may be represented by a new region 
adjacency graph RAG2, whose nodes are obtained by contracting each tree 
of the forest MSFi in the graph RAGi. Repeating the same treatment to 
the graph RAG2 produces the next level of the hierarchy, where each tree 
of the minimum spanning forest MSF2 has been obtained by merging 
several trees of the MSFi. The process is then repeated until a graph is 
created with only one regional minimum. 

' r 1 Ill I mill I II I II II I 



Fernand MeyerCentre de Morphologie MaXhe Watersheds, waterfalls, on edge or node weigl 2012 February 29 187 / 201 



Construction of the level 1 of the hierarchy 



We start with an arbitrary node or edge weighted and connected graph 
G. As explained earlier, we extract a graph flooding graph G' . For a 
steepness k, we prune G' and get J,'' G/. The scissor operator x produces 
a minimum spanning forest Fi = x Gf, spanning the finest watershed 
partition, the lowest level of the hierarchy. 

The graph representing the second level of the hierarchy is obtained by 
contracting all edges of the forest Fi ; each tree becomes one node. The 
nodes are connected by edges of the graph G. The result of this 
contraction G[ = k{G, Fi). is again a connected graph, to which the same 
treatment as previously can be applied. 
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Construction of the level 2 of the hierarchy 



Only the nodes of the graph G[ are weighted : G[ = (ei,o). The nodes 
will be weighted with £„£ and we get the graph {ei,£neSi), which becomes 
a flooding graph of steepness k in J,'' (ei, Ene^i)- A final scissor operator 
creates a forest F2 spanning nodes of G[. F2 represents the level 2 of the 
hierarchy, and the nodes of G spanned by each of its trees constitutes a 
catchment basin of level 2. 

The contraction k(G(, Fi) produces again a connected graph = (ei,o) 
to which the same treatment may be applied. 

This process is repeated until the graph k{G'^, Fm) contains a unique 
regional minimum of edges. 



I' 1 II I llllll I II I II II I 



Fernand MeyerCentre de Morphologie MaXhe Watersheds, waterfalls, on edge or node weigl 2012 February 29 189 / 201 



Emergence of a minimum spanning tree 



Each new minimum spanning forest makes use of new edges of the 
initial graph. The union of all MSF constructed up to the iteration m, 
U F^t, is a minimum spanning forest of the initial graph G, which 

k<m 

converges to a MST of G as m increases. Hence the union of all edges 
present in the series F^ is a minimum spanning tree of the graph G. 
The particular minimum spanning tree which emerges depends upon the 
depth k of the pruning operator J,'*, and of the particular choices among 
alternative solutions made by the pruning operators x used at each step. 
We call the operator extracting the minimum spanning tree }i{G). 
The figure in the next slide presents a RAG, its flooding graph in the first 
row, and in the second row the MSFs |J for m = 1,2, 3. The last one 

k<m 

being the MST. 
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Emergence of a minimum spanning tree 




Emergence of a minimum spanning tree 



If G and G' are two flooding trees, G' being a partial tree of G included in 
G, then the pruning operators x have less choices for pruning Gf as for 
pruning G. For this reason the family of MST derived for G' is included in 
the family of MST derived for G: {li{G')} C {f<(G)} . 
In particular if one considers the decreasing sequence of minimum 
spanning forests x ^-^ G, one obtains decreasing families of MST : for 
/c < /, we have {/^(i' G)} C G)} . The choices are more and more 

constrained. One gets minimum spanning trees which are steeper and 
steeper for increasing k in J,'' G. 
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Discussion 



All MSTs of an edge weighted graph share a fundamental property : 
between any two nodes of the graph, there exists a unique path in the 
MST which links them, and this path has a minimal flooding weight. 
Thanks to this property, replacing the graph G by its MST T permits to 
construct the catchment basins linked to the lexdisti . This procedure is 
fast but has no control on the quality of the result if one chooses an MST 
at random (the one which is produced by the preferred MST extraction 
algorithm). Fast, as there are as the number of edges is strongly reduced 
(in a tree : number of edges = number of nodes -1). Limiting, as one has 
no control on the steepness of the watershed which is ultimately extracted. 
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Discussion 



In order to get high quality partitions and segmentations, one has to 
carefully chose the MST. The preceding slide has shown that the MSTs 
form a decreasing family of trees as their steepness increases. Higher 
quality segmentations will be obtained for a higher steepnes. In order to 
improve things, one may prune the graph G with the operator before 
extracting one of its MST, yielding a MST of k-steepness. 
One has then to consider a trade-off between the required speed and the 
quality of results one targets. This obviously also depends on the type of 
graphs. Some graphs contain only very few MST ; in the extreme case, a 
graph where all edges have distinct weights, contains a unique MST. 
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Conclusion 
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Outcome 



Starting with the flooding adjunction, we have introduced the flooding 
graphs, for which node and edge weights may be deduced one from the 
other. 

Each node weighted or edge weighted graph may be transformed in a 
flooding graph, showing that there is no superiority in using one or the 
other, both being equivalent. 

We have introduced pruning operators J,'* and ^^^^ which extract 
subgraphs of increasing steepness. For an increasing steepness, the number 
of never ascending paths becomes smaller and smaller. This reduces the 
watershed zone, where catchment basins overlap. 
The scissor operator x, associating to each node outside the regional 
minima one and only one edge choses a particular watershed partition. 
Again, with an increasing steepness, the number of equivalent solutions 
becomes smaller. Ultimaterly, for natural image, an infinite steepness leads 
to a unique solution, as it is not likely that two absolutely identical non 
ascending paths of infinite steepness connect a node with two distinct 
minima. Fl 1IIIIIIIII II I II II I 
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Outcome 



Finally, we have shown that the NAP paths remaining after the pruning J,'' 

or ^(''^ are identical with geodesies of lexicographic distances. 

We have shown how the path algebra may be adapted to lexicographic 

distances. 

We have presented the Moore-Dijkstra algorithm for constructing skeletons 
by zone of influence for lexicographic distances. 

For the depth 1, the lexicographic distance becomes the flooding distance 
and the algorithms become classical algorithms for constructing minimum 
spanning forests. 

For the depth 2, it is equivalent with the topographic distance. 
The beneficial effect of hierarchical queue algorithms is highlighted, as it 
permits to correctly divide plateaus among neighboring catchment basins. 
We also presented the core expanding algorithms which are particularly 
efficient for lexicographic distances. 
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Outcome 



The waterfall hierarchy is obtained by contracting the trees of a first 
watershed construction and the new graph submitted to a novel watershed 
construction. The process is iterated until only one region remains. 
The union of the edges of all forests produces constitute a minimum 
spanning tree of the initial graph. 

MST and MSF defined before only for edge weighted graphs may be 
extended to node weighted graphs. 

MST and MSF may be ordered into nested classes according to their 
steepness. 

The classical order is thus inverted : classically a MST is extracted from 
the graph as a first step, and a MSF derived from it, by cutting a number 
of its edges. This way of doing leaves not much control as the quality of 
the result depends upon the choice of the MST. 

As MST and MSF may be ordered into nested classes according to their 
steepness. 
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Outcome 



The algebraic approach to the watershed presented here sets the stage in 
which a number of earlier definitions and algorithms may be reinterpreted : 

• the waterfall algorithm P| 

• the watershed line defined a zone of influence of the minima for 
various distances, in particular the flooding distance and the 
topographic distance [ll],[23]. In the case of the flooding distance 
applied on an edge weighted graph, one finds the graph cuts [5], and 
the algorithm for constructing a waterfall hierarchy described in [16] 

• the role of the hierarchical queues for a correct division of plateaus 
between catchment basins if one uses myopic distances |l7j 

We also explored the place of the choice for constructing a watershed 
partition. The choice becomes more and more restricted as one considers 
lexicographic distances with increasing depth. 
The waterfall hierarchy transposed to MST provides a hierarchical 
decomposition of the MST. ^ ^^^^^^ ^ ^ 
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Outcome 



Old algorithms may be reinterpreted, in particular those relying on pruning 
and labeling graphs [7], [15], [24], [25], 

The classical algorithm for constructing watersheds derived from the 
hierarchical queues is a particular case of core expanding algorithm. 
In order to reduce the number of equivalent watershed partitions one is 
able to extract from a flooding graph, one may use a mixture of pruning 
^(''^ up to a given depth, and on the resulting graph apply a myopic 
algorithm, with a lexicographic depth of 1 or 2. 
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The watershed transform introduced by S.Beucher and C. Lantuejoul [4] is 
one of the major image segmentation tools, used in the community of 
mathematical morphology and beyond. If the watershed is a successful 
concept, there is another side of the coin: a number of definitions and 
algorithms coexist, claiming to construct a wartershed line or catchment 
basins, although they obviously are not equivalent. We have presented 
how the idea was conceptualized and implemented as algorithms or 
hardware solutions in a brief note :" The watershed concept and its use in 
segmentation : a brief history" (arXiv:1202.0216vl), which contains an 
extensive bibliography. Here we give a more restricted list of references. 
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